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Abstract 



Laplace transform QCD sum rules for two-point functions related to the strangeness- 
changing scalar and pseudoscalar Green's functions ip{Q'^) and ip^iQ'^), are used to 
determine the subtraction constants ■0(0) s-^id ^^5(0)) which fix the ratio Rgu = <f|>- 
Our results are V'(O) = -(1.06±0.21) x IQ-^ GeV^ 05(0) = (3.35±0.25) x lO^^ QeV^, 
and Rsu = = 0.5 it 0.1. This implies corrections to kaon-PCAC at the level of 

50%, which although large, are not inconsistent with the size of the corrections to 
Goldberger-Treiman relations in SU{3) ^ SU{3). 
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The quark vacuum condensate ratio Rsu = ^'^u> ' < uu > ~< dd >, is an important 
quantity which measures the size of flavour SU(3) symmetry breaking in the QCD vacuum 
[|1|]. It also enters in a variety of QCD sum rules 0, particularly in those used to determine 
baryon masses 0. In addition, it provides a measure of the corrections to kaon PCAC 
(Partial Conservation of the Axial- Vector Current), provided the quark masses are known 
and neglecting corrections to pion PCAC In fact, let us consider the two-point 

function involving the flavour-changing axial-vector current divergences 



^5(g') =ijd^x e^'i^ (0|T(9M^(x) a'^4(0))|0; 

where flavour indices have been omitted for simplicity, and 

d^'A^{xX = {rrii + mj) : ijj{x) ^75 ^x) : , (2) 

and (i,j) are flavour indices. A well known Ward identity fixes the value of ip5{<f) at zero 
momentum 0,0], viz. 

^5(0)1- = -(mi + rrij) l^^^^ + , (3) 

which is a renormalization group invariant quantity. Saturation of Eq.(l) with the lowest 
lying pseudoscalar meson then leads to the Gell-Mann, Oakes, Renner relations 

- (m„ + mrf) (mm + 3c?) = 2fl^l , (4) 

for (i,j) = (u,d), and 

- {rris rriy) {ss + uu) = 2f'^^\ , (5) 

where now (i,j) = (u,s), = 92.4 ± 0.4 MeV, and fx = 113.0 ± 1.5 MeV §. Considering 
Eq.(3) for up- and strange-quark flavours and then for up- and down-quark flavours and 
taking the ratio leads to 

where < uu > ~< dd > is a very good approximation in this case 0. In principle, this 
relation could be used to determine Rgu, if the ratio of the subtraction constants '?/'5(0), i.e. 



Raa^ is known independently. Using current values of the quark masses in Eq.(6) leads 
to 

< e c > 

R,^ = — ~ 0.15 Raa-1. (7) 

< uu> 

Since Raa is expected to be of order (9(10) (from using PCAC), this procedure would lead 
to very large uncertainties, unless V'slO) were to be determined with extreme accuracy. An 
alternative method not affected by this problem was first proposed in [0, and it is based 
on examining ^/'5(0) together with '?/'(0), where 



V;(g2) = i j (fx e'"" {Q\T{d^'V^{x) 9"K^(0))|0 



with 

d^V^{x)\l = {rrij - rrii) : tjjjix) i il)i{x) : , (9) 

being the flavour-changing vector current divergence. The Ward identity analogous to 
Eq.(3) is now 

^(0)1 = -(m, - m,) (^^.^, - ^,^,) . (10) 

Defining the ratio 
it follows that 

fl,„.M^i±|i:i. (12, 

{uu) 1 - RvA 

where the reasonable approximation (m^ — m^) ^ (m^ + m„) has been made. This pro- 
cedure can provide a more accurate determination of Rsu, as discussed in 0. Since then, 
considerable progress has been made on the QCD calculation of the two-point functions 
Eqs.(l) and (8), viz. the perturbative expressions are now known up to four loops [§, 
quark-mass corrections have been calculated up to order (9(m^), and the issue of quark- 
mass logarithmic singularities has been satisfactorily resolved [@|-llO|. On the other hand, 
there is now better experimental information on the hadronic spectral functions in both 



the scalar |11| and the pseudoscalar channels pf. In view of this, we discuss here a new 



determination of the subtraction constants '?/'(5)(0) which are then used to determine the 



3 



ratio Rsu through Eq.(12). 



In order to relate the constants ip[5){0) to the two-point functions 4'(5){q'^), Eqs.(l) and (8), 
one defines the auxihary functions 

g2D(5)(g2)=^(^)(^2)_^(^^(0). (13) 

The derivatives of these functions, i.e. 

= ^^(?') , (14) 

<Pi<i') = ^^5(g^) , (15) 

have the same imaginary part as il'{5){Q'^)y ^ind satisfy the following dispersion relations 

e(Q') 1 _ 1 rds Irmlj(^5){s) ^^^^ 



?') j Trio s (s + g2)2 ' 
where = —q^ > 0. After Laplace improvement, these dispersion relations become 

e(Mi) 1 _ 1 rds 



0(M2 



e-^/^-^Wm^(5)(s) . (17) 

L) \ IT JO s 



According to quark-hadron duality, the left hand side of Eq.(17) is computed in QCD 
using the Operator Product Expansion to organize the perturbative series and the power 
corrections, and the right hand side is saturated with the experimental data. Using the 
available QCD expressions for il'i5){(l'^) IH'IUHI? performing the derivatives, Laplace trans- 



forming, and after Renormalization Group improvement, we obtain the following QCD 
representations for the left hand sides of Eq.(17) for three quark fiavours 
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where '^{z) is the di-gamma function, primes stand for its derivatives, and C,{z) is Riemann's 
zeta function. The four-quark vacuum condensate term of dimension-six has been omitted 
above as it makes a neghgible contribution, and the up-quark mass has been neglected in 
comparison with nis- The four-loop expressions for the running coupling and quark mass 



in the MS renormalization scheme, and for three flavours, are 
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where L = ln{Q'^/AQ^,j^), LL = InL, and 

281198 



/?4 



4608 



fc(3), 



(21) 



with 74 = 88.5258 |T2[, and mj is the invariant quark-mass. The terms of order (^(x*) 
above are known up to a constant not determined by the renormahzation group. This 



constant can be estimated e.g. using Fade approximants |T3[. However, we have checked 
that our final resuhs are essentially insensitive to terms of this order. 



Before we proceed to discuss the hadronic parametrization of the scalar and pseudoscalar 
spectral functions, it should be stressed that these spectral functions are also used to de- 
termine the strange-quark mass from QCD sum rules for ip5{q'^) [0 and iplq"^) p|]-[p^,|r5|. 
The fact that is one of the most important parameters in Eq.(18), and given the wide 
range of values obtained from QCD sum rules [0, it is imperative to achieve consistency 



between the results for from the scalar and the pseudoscalar channels. Fully consistent 



results are obtained as explained below. Another very important parameter in Eq.(18) 
is Aqcd- The value extracted from a variety of experimental data has been increasing 
steadily over the years, from Aqcd — 100 — 200 MeV at the birth of the QCD sum rules 
0, to something as high as kqcoiNf = 3) ~ 300 — 450 MeV lately 0. Such high values 
make radiative corrections to Green's functions the overwhelming terms in the Operator 
Product Expansion, thus reducing the numerical importance of the power corrections as 



parametrized by quark and gluon vacuum condensates. In fact, it has been argued in |]I7 



that if AqcniNf = 3) exceeds ~ 330 MeV the QCD sum rule program may break down, 
and it becomes extremely difficult to extract numerical values for the condensates from the 
data (see also JT^). As for the other QCD parameters in Eq.(18), the gluon condensate 
has the value [|19[ (^G^) ^ 0.01 - 0.04 GeV^ The light-quark condensate (at 1 GeV) 
is (uu) ~ —10^^ GeV^. For the strange-quark condensate, since it is the object of this 
determination, an iteration procedure of fast convergence has been adopted, starting from 
(ss) ~ {uu). 



Turning to the hadronic sector, the spectral function in the pseudoscalar channel is parametrized 
as in 113, viz. 



-l7m/j,is)\HAD = 2f'^M],6 is - 

71 ^ 
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where the threshold behaviour of the resonant piece, determined using chiral perturbation 
theory, is given by 
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(24) 



7 



where 

, _ {Ml: - Mj-f + Ml-.Tj. 

is the contribution from the resonant sub-channel K*{892)—7t. The parameterA above con- 
trols the importance of the second radial excitation (K(1830)) relative to the first (K(1460)) 
(A ^ 1), with both resonances being parametrized by Breit-Wigner forms BWi^2{s)- Due 
to the approximation rriu = 0, the pion mass has been neglected. There is another resonant 
sub-channel, the p{770)K, which turns out to be numerically negligible [14|. As usual, this 



resonant hadronic parametrization is used up to some threshold energy, sqa, the continuum 
threshold, after which the spectral function is assumed to be given by perturbative QCD. 
In the case of the scalar channel, there is experimental data on Kit phase shifts []TT| which 
can be used to reconstruct the spectral function given by 



1 3 a/(s - s+) (s - S-) „ 

-Im ^ [s) = -f^ tl^ L \d{s)f , (26) 

where s± = [fix ± A^tt)^, and d{s) is the scalar form factor. We have used the method of 



| 15| , based on the Omnes representation, to relate d{s) to the experimental phase shifts. In 
this way, the results for from the pseudoscalar channel are fully compatible with those 
from the scalar channel, giving an invariant strange-quark mass 

= 140 ± 10 MeV , (27) 

for Aqcd in the range Aqcd — 300 — 350 MeV. As mentioned earlier, higher values of Aqcd 
may invalidate the QCD sum rule program in general, and they do so in this particular 
application, as they lead to serious instabilities in the results for '?/;(5)(0). Therefore, we 
restrict Aqcd to the above range. In Figs. 1 and 2 we show typical results for ip^iO) 
and V^(0), respectively, obtained using the values Aqcd = 300 MeV, = 145 MeV, 
(^G^^ = 0.01 GeV^, and the asymptotic freedom thresholds sqa = 6 GeV^, and sov = 
4 GeV^. These results lead to the ratio Rsu shown in Fig. 3. As it may be appreciated from 
these figures, the results are nicely insensitive to changes in the Laplace variable M|; they 
are also reasonably stable against changes in the continuum thersholds, sqv,a- Increasing 
Aqcd leads to smaller values of m^, and higher values of the ratio Rgu- After allowing for 
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variations in all the relevant parameters, within the ranges indicated above, we obtain 

^5(0) = (3.35 ± 0.25) X IQ-^ GeV^ , (28) 

^(0) = -(1.06 ± 0.21) X 10-3 GeV^ , (29) 

R,^ = = 0.5 ± 0.1 . (30) 

< uu > 

A comparison of these results with previous determinations based on similar methods 0, 
pO| is not very enlightening, as most of the old analyses were done at the two-loop level in 
perturbative QCD, they were affected by logarithmic quark-mass singularities, and did not 
make use of the latest experimental data to reconstruct the hadronic spectral functions. 
However, these problems do not affect determinations based on different methods. For 
instance, from QCD sum rules for baryon masses |^ Rsu = 0.6 ±0.1 has been obtained. 



while an estimate based on non-perturbative quark-mass generation gives Rsu — 0.7 
(with no error estimates). Our result, Eq.(30) is consistent with these values. Furthermore, 
Eq.(29) leads to a correction to kaon-PCAC, 6k, defined as 

^,{0) = 2flixU^-5K), (31) 

where using Eq.(29) one finds 

5k ^ 0.5 . (32) 

This result points to a rather large correction to kaon-PCAC, although it is not inconsistent 
with the expected size of the corrections to Goldberger-Treiman relations in SU {3)^SU (3) 
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Figure Captions 

Figure 1. The subtraction constant ip^i^) a function of the Laplace variable M|, for 
Aqcd = 300 MeV, = 145 MeV, (^C^) = 0.01 GeV^ and the continuum thresholds 
SoA = 6 GeV^ and Sov = 4 GeV^. 

Figure 2. The subtraction constant ip{0) as a function of the Laplace variable M|, for 
Aqc^ = 300 MeV, m, = 145 MeV, (^G^) = 0.01 GeV^ and the continuum thresholds 
sqa = 6 GeV^, and sqv = 4 GeV^. 

Figure 3. The ratio Rsu = ^f^^ as a function of the Laplace variable M|, for Aqcd = 
300 MeV, nis = 145 MeV, (^G^) = 0.01 GeV^ and the continuum thresholds sqa = 
6 GeV^ and sov = 4 GeV^ 
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